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Abstract 

We study four-dimensional simplicial gravity through numerical simulation 
with special attention to the existence of singular vertices, in the strong cou- 
pling phase, that are shared by abnormally large numbers of four-simplices. 
The second order phase transition from the strong coupling phase into the 
weak coupling phase could be understood as the disappearance of the singu- 
lar vertices. We also change the topology of the universe from the sphere to 
the torus. 
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1 Introduction 



One of the most exciting challenges in the theoretical physics is to understand gravi- 
tational interaction in the context of the quantum theory. The problem we encounter 
when we try to formulate quantum gravity within ordinary field theory in four di- 
mensions is that we cannot renormalize it perturbatively. If we use lattice regular- 
ization, which enables a nonperturbative study, general coordinate invariance is not 
manifest and whether it is restored in the continuum limit is a crucial problem. One 
possibility of lattice regularization of quantum gravity is dynamical triangulation, 
which is believed to restore general coordinate invariance in the continuum limit. 
It has been solved exactly in two dimensions and its continuum limit is shown to 
reproduce Liouville theory, in which general coordinate invariance has been treated 
carefully. 

Although four-dimensional dynamical triangulation seems to be difficult to solve 
analytically, there is no potential barrier in studying it through numerical simula- 
tion. Employing the Einstein-Hilbert action as the lattice action and sweeping the 
gravitational constant, it has been discovered that the system undergoes a second 
order phase transition g, which suggests the possibility of taking a continuum 
limit. One of the main purpose of this paper is to try to clarifying the physical 
meaning of this phase transition. 

When we regularize four-dimensional quantum gravity with dynamical triangu- 
lation the integration over the metric is replaced with the random summation over 
all possible four-dimensional simplicial manifolds. 

Although we can modify the lattice action expecting universality, it is natural 
to start with the Euclidean Einstein-Hilbert action 

s = J dVs (a - ^r) (i) 

as a first trial, where A is the cosmological constant and G is the gravitational 
constant. Let us denote the number of i-simplices in a simplicial manifold by iVj. 
One can easily find that, for a simplicial manifold, 

fd*xy/g = cN, (2) 

J d 4 x^gR = 2itN 2 - 10a-/V 4 , (3) 

where c is the volume of each four-simplex and a is the angle between two faces of 
a four-simplex, which is equal to arccos y^j. Therefore the Einstein-Hilbert action 
(HD can be expressed in terms of lattice variables as 

Slat = K 4 iY 4 - k 2 N 2 , (4) 



where k 4 and k 2 are related to A and G through, 

10a 2tt 



2 The Vertex Order Concentration 



We consider an ensemble with a fixed N4 and with spherical topology. There are well 
established methods for generating such an ensemble through numerical simulations, 
and the technical details of our simulation shall be given elsewhere . Our code is 
written for arbitrary dimension following Ref. Q. 

Let us turn to the results of our simulation. We first look at the second order 
phase transition, which can be seen through thermodynamic quantities such as the 
average curvature per unit volume : 

i?av = ^ (Rtot = I d'x^R) (6) 



i 



iV 4 

Fig. |l| shows our results for (R av ) at various fi^'s 



2tt— 4 - 10a. (7) 




Figure 1: The average curvature per volume is plotted against k 2 for N4 = 
8000, 32000. 



In contrast to the three-dimensional case [||, no hysteresis has been observed. 
Also, one sees that the size dependence of the data changes abruptly at k 2 = 1.2 ~ 
1.3. On the right there is little size dependence, whereas on the left, the curve 
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goes lower and lower as we increase the system size. The derivative of the average 
curvature gives the susceptibility 



XR = 27T 



d(R ay ) « t > - (R tot y 



8k 2 iV 4 



which represents the fluctuation of the total curvature. As is expected from Fig. [l|, 
the susceptibility has a peak around k 2 = 1.2 ~ 1.3, which grows higher as the sys- 
tem size is increased. This implies that the correlation length of the local curvature 
diverges at the critical point ||, where we may hope to take a continuum limit. 
Since k 2 corresponds to the inverse of the gravitational constant, as is seen from (||), 
we call the large k 2 phase as the weak coupling phase and the small k 2 phase as the 
strong coupling phase . 

Although this phase transition has been observed by many authors]!], the 
physical origin of this transition might be not understood clearly. 

To clarify it, we measure the vertex order distribution as follows. Let the vertex 
order o(v ) be the number of four-simplices sharing the vertex v. Then the vertex 
order distribution p can be defined as 

p(n) = -^r(Y, 5 o(v),n)- (9) 

This quantity is measured every 100 sweeps and averaged over 100 configurations. 
In order to reduce the fluctuations of the distribution, we smear the data over bins 
of size 10. 

We first note that the average vertex order per one vertex, o(v), can be given as 

o(v) = — > o[v) = = . 10 

1 } N V N N 2 -2N 4 + 4 1 ; 

In the second equality, we used the relation J2 V °i v ) = SN^, which comes from the 
fact that each four-simplex has five vertices. 

From Eq. ([|) and Eq. (|1(]), when we move from the strong coupling phase (^2 ~ 
0.0) to the weak coupling phase (k 2 ~ 2.0) with fixed system size, A^ 2 increases and 
thus o(v) goes to a smaller value. 

In Fig. H we show the vertex order distribution p(n) for k 2 = 0.0, 1.267 (near 
the critical point) and 2.0 with N 4 = 32, 000. For k 2 = 0.0, one finds that there 
is an isolated peak of very large vertex order; as large as one third of the total 
four-simplices. 

As k 2 is increased from 0.0 the position of the peak shifts to left in accordance 
with the decrease of the o(v) and around the critical point (k 2 ~ 1.2) the peak 
is absorbed into the continuum part of the distribution. For k 2 = 2.0, on the 
other hand, the distribution damps quite rapidly for large vertex order and we have 
confirmed that it remains almost unchanged when we increase the system size. 



4 




Figure 2: The vertex order distribution for k 2 = 0, 1.267, 2.0 for iV 4 = 32000. 



In Fig. [|we show the size dependence of the vertex order distribution for k 2 = 0.0. 
One finds that the very large vertex order grows linearly as one increases iV 4 , and 
thus this concentration might be relevant even in the thermodynamic limit. 

We call this peculiar phenomenon as the vertex order concentration (VOC). We 
have also confirmed that the peak consists of two vertices^, and we found no link 
order concentration; there is no singular link shared by conspicuously large numbers 
of four-simplices. 

1 For D-dimensional dynamical triangulation(3 < D < 7), it seems that the peak consists of 
D — 2 vertices. 
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Figure 3: The vertex order distribution for k 2 = varying the system size iV 4 from 
8k to 128k. 



3 The thermalization check 

One may suspect that our Monte Carlo simulation is trapped by a local minimum 
(meta-stable state) and sweeps over only a limited part of the whole configuration 
space. To check that this suspicion is not the case, we prepare three types of 
configurations for initial configurations from which we start simulations. 

The first one consists of six four-simplices which are the surface of a five-simplex; 
we call this the hot start configuration. The second one is the cold start, for which we 
prepare approximately flat configuration, using the surface of the five-dimensional 
rectangular complex (box). The system size of the cold start configuration can be 
adjusted near the target number N® and it has no VOC[Q]. The third configuration, 
four-VOC configuration, has four singular vertices and is made in the following way. 
We prepare two configurations of system size iV 4 /2 by performing sufficiently many 
sweeps in strong coupling phase. Each of the thermalized configurations has two 
singular vertices. Then we identify one four-simplex of one configuration with one 
four-simplex of the other, so that the resulting manifold is the four-sphere of system 
size iV 4 — 1 with four singular vertices. So this configuration has twice as many 
singular vertices as the configuration obtained in strong coupling region. 
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After performing more than 10,000 Monte Carlo sweeps [] we found no depen- 
dence on the three types of initial configurations for any of our observables. 

For example, we show the history of the vertex order below. Let us label the 
vertices of the configuration as Uj, so that o(i>j) > o(t>i + i) is satisfied for any i = 
1, 2 • ■ ■ iVo- We show the o(t>i) (the largest vertex order) as a function of the number 
of Monte Carlo sweep from each initial configuration together with 0(^2) and 0(^3) 
for K2 = 0.0 in Fig. |j. We only show the results of the simulation for the cold and hot 
starts for the legibility. For each types of initial configuration (thick curves or thin 
curves), one can see that there are two large vertex orders o{vi) and 0(^2) around 
25, 000 and a large gap between o(vs) and the above two. Below the curve of 0(^3), 
the curves of o(vi),i > 3 run closely without significant gaps. One can easily see 
that the two curves of each o(t> j), i = 1, 2, 3 are almost identical after 10,000 sweeps. 

N 4 =7488 k 2 =0.0 

3000 1 




1 • 1 • 1 

10000 20000 

sweep 

Figure 4: The history of the vertex order for n 2 = for the hot and cold initial 
configurations. 

Considering that the vertex order distributions of the three types of configurations 
are quite different from each other, the above fact strongly suggests that our result, 
the existence of two singular vertices in strong coupling region, does not come from 
the insufficiency of thermalization but reflects true properties of the path-integral 
measure of the dynamical triangulation. 

2 we define one sweep as N% times accepted updates. 
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4 The torus topology 



We test a possibility of the phenomenon of the VOC being related to the constraint 
of the topology of the manifold, i. e the manifold must be S A . We change the manifold 
from four-dimensional sphere (S" 4 ) to the torus ((S" 1 ) 4 ). The method of making the 
(S* 1 ) 4 triangulation is as follows. 

A) Prepare the four- dimensional rectangular complex (four-box) and identify each 
pair of parallel boundary. Draw one line Id between two vertices which is 
diagonal to each other. 

B) Divide the four-boxes into 4! four-simplices. The dividing edges on boundaries 
of the four-boxes must be projections of the diagonal line of Id for consistent 
construction of the torus triangulation. 

For simplicity Fig. ||] describes the three-dimensional case. 
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Figure 5: The construction of the torus triangulation. 



The resulting triangulation has (S 1 ) 4 topology, from which we start the Monte 
Carlo simulation. 

Fig. |6| shows the comparison of the vertex order distribution p(n) between the 
sphere and the torus case for A4 = 128, 000. The p(n) of the torus of the size 
A4 ^ 100, 000 is sensitive to whether the manifold is the torus or the sphere, while 
for the larger manifold A 4 ~ 128, 000 the distribution is almost identical between 
the two topologies (Fig. |5|). 
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Figure 6: The vertex order distribution of the sphere and the torus. A/4 = 
128, 000, k 2 = 0.0. 

The number of the singular vertices in torus manifold is one for iV 4 £ 100, 000 
and two for A/4 ~ 128, 000. From these results we might say that the finite size 
effect of torus is severer than the sphere case. And the vertex order distribution 
p(n) seems to be almost insensitive to the topology difference, but of course more 
high statistic simulations of more variation of topology (e.g. (S 2 ) 2 , D A etc.) are 
necessary for the definite conclusions. 

5 Discussion 

To summarize, we found two singular vertices in the strong coupling region, which 
disappear in weak coupling region. Although the thermalization is checked, the 
existence of this VOC is still very strange. The VOC in k 2 = case, for which we 
observe purely the measure of path integral without any weight from the action, 
implies that the number of such VOC configuration is much larger than the number 
of the smooth (non VOC) configuration in dynamical triangulation. Such a singular 
behavior might be an obstacle to the continuum limit. Considering universality in 
quantum gravity, as well as in ordinary field theories, we think that a sound second 
order phase transition without VOC, where we can take a sensible continuum limit, 
should be searched by modifying the lattice action]?], |2|, || to suppress the VOC. 
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